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Abstract 

We give a simple expression in linear and quadratic Dalitz-plot slopes 
which does not depend on the charge combination of the 3ir state (K ± — > 
7r ± 7r + 7r _ or 7r ± 7r°7r c) and K® — > tt + tt^tt° or 7r 7r°7r°), if all phases between 
final states are negligible. After investigating the influence of radiative 
corrections, it is shown how new measurements, especially of quadratic 
slopes in the ir^^iT channel, could help to test theoretical predictions 
more stringently. 







Introduction 



At least until the advent of decisive data on CP-violation from the B sector, 
the K sector remains our main source of experimental information on this old 
and important question. Further progress depends, however, not only on very 
high statistics experiments, but also on well tested effective theoretical models 
for handling long-range strong interactions in a quantitative way. 

As in the 3-generation quark model CP-violation is thought to be due to a 
phase factor in the Cabbibo-Kobayashi-Maskawa-matrix, any measurable first 
order effect (,, direct" CP-violation) is expected to disappear if certain phase 
shifts from (strong) final state interactions (f.s.i.) between different components 
of transition amplitudes are vanishing. In an attempt to estimate the size of 
CP-violation effects in 3ir decays |IJ it was found necessary not only to 
fit a number of coefficients determining the weak effective Lagrangian from 
experimental data, but also to rely on loop calculations in deriving the f.s.i. 
phases. A further problem is that the strong interaction effects in 3ir 
decays cannot be handled (like in K— * 2tt) by just multiplying final (strong) 
eigenstates with constant phase factors; imaginary parts from pion loops depend 
on dynamical variables, and the real parts are modified in turn by various loop 
contributions, among them loops with inner K lines. Therefore the denotation 
as f.s.i. phases is not quite correct. It would be desirable to get more direct 
experimental information on these phases from existing experimental data on 
K decays, but this turns out to be very difficult, as we will explain below. 

A well documented joint fit of isotopic K— ■> 2ir, 3tt amplitudes has been 
given by Devlin and Dickey and more recent work is published by Kambor 
et al.||. Both groups of authors found it necessary to make assumptions on the 
f.s.i. phases in 3tt states, namely that they are zero, resp. < 15°. While general 
unitarity arguments lead us to agree with these expectations concerning the 
constant parts of the amplitudes, the situation may be more complex for the 
coefficients (bu, cu, du, see below) of the non constant parts, which vanish at 
the center of the Dalitz plot. 

Instead of repeating once more these fits, varying merely the data sets used 
and the parameter sets to be fitted, it may be more useful to investigate by 
simpler means which data are particularly important for our question, and 
to check their internal consistency and the appropriateness of the radiative 
corrections applied before the data can be used to fit isospin amplitudes. Their 
effect is found to be very important for our investigation, and the standard 
procedure for handling them (see Q) may not be sufficient. 

In the following sections, we first repeat the definitions and notations (sec- 
tion 1), discuss the internal consistency of slope parameters (section 2), in- 
troduce radiative corrections and demonstrate a method to recalculate slope 
parameters without reference to raw data from older experiments (chapter 3), 
and finally discuss the significance of future experiments (chapter 4). 
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Definitions and Notation 



K+ - 


-> 7T + 7T + 7r : 


A+ 


K+ - 


-> 7r°7r°7r+ : 


A+' 


Kl- 


-> 7r + 7r _ 7r° : 


A 


Kl- 


-> vr 7r 7r° : 


A ' 



In order to fix the notation, we give a short account of relevant definitions for 
3-7T decays connected to isospin relations. 
The kinematic (Dalitz plot) variables used are 

y _ (gg ~ gi) v _ («3 - ap) , 

2 ' 2 ' 

«i = (PK ~ Pi) 2 (i = 1, 2, 3) , s = -(si + s 2 + s 3 ) • 

We further use the three (not independent) variables rj = s« — so for expansions 
about the center of the Dalitz plot. We will consider the following transition 
amplitudes (K — > ttitt^s): 

: 2~2 A tj(n,r 2 ,r 3 ) , 
--H A fj{n,r 2 ,rz) , 



Factors (n!) -1 / 2 for identical pions are not included. Each A is a sum of isospin 
contributions with / = final state isospin, contributing to the given channel, and 
J = two times the isospin change AI. The relevant states have I = 1 or I = 2 
for isospin changes AI = 5, |- The (37r)j=o-state is totally antisymmetric and 
therefore has a totally antisymmetric momentum eigenfunction, that means a 
form factor / ~ (r\ — r 2 )(r% — r 3 )(r 2 — r 3 ), leading to contributions of at least 
3 rd order in rj, which we will neglect. 1 = 2 final states are not present in 
decays, if we assume CP invariance (CP = (— 1) for 37r S-states). 

Due to the Wigner-Eckhart-theorem, we have to introduce only 3 inde- 
pendent form factors /n,/i3,/23, whereby further restrictions are derived by 
expanding the amplitudes up to 2 nd order and separating different symmetry 
classes in r\ , r 2 , r 3 . 

Defining / (1) (r 1; r 2 , r 3 ) = \[f{n, r 2 , r 3 ) + f(n, r 3 , r 2 )], etc., the expansion 
up to 2 nd order in ri )2 , 3 can be written: f®(ri, r 2 , r 3 ) a + 6rj + c(rf + 
r l + r l) + rf(2r? — r ] — r k) (other linear and quadratic terms can be included 
by redefinition of b,c,d, since ri + r 2 + r 3 = 0). After constructing properly 
symmetrized I = 1,2 final states from three pion (Ij = 1) states (sec || for 
details), one finds 

A+ _ f (l) , f (2) 
A +' = 

= (J = 1,3), 

4+ _ j+'_n/(3) f (l) f (2) 
^23 — ^23 — z /23 J23 J23 
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with ^ 

-gu - V2 ^ .J u ~ {-2} fu I 1 ) 

V2'~ l ~2'2'~ , ~2 I > / 

for J = 1, 3 (the factor y/2 follows from K\ = l/y/2(l - CP)K°). 
Expanding ffj in the above manner and using 

r 3 ~y, r 2 + r 2 + r 2 ~Y 2 + X 2 /3 , 2r§ - r 2 - r 2 2 ~ Y 2 - X 2 /3 , (2) 
the four amplitudes can be written as 

^ = + + c "( y2 + X V3) + a/j(^ 2 - X 2 /Z)\ (3) 

with the coefficients given in table l.(we write, with slight redefinition, again 
E«u = a etc.) 

The (irrelevant) relative (— )sign between K° and X + amplitudes has been 
chosen in accordance with 0. We do not distinguish at this stage between 
charge conjugated channels and take in the following also averaged exper- 
imental data for them. 

Due to strong interactions, that means rescattering between initial and/or 
final states, the fu become complex functions of r$; therefore we have to consider 
relative phases between all coefficients. Usually an is chosen to be real and 
positive. 

Consistency of Slope Parameters 

From the usual definition of the measurable slope parameters g, h, k (in each 
channel) : 

\A\ 2 ~l + gY + hY 2 + kX 2 

one gets 

_ Re(ab*) u _ | b | 2 +2Re[a{c* + d*)] , _ 2 Re[a{c* - d*)\ 
9 — 2 j , h — j , k — — j . 

\ a \ z | a | z 3 | ° | 

Since, as mentioned in the introduction, the strong phases for the constant terms 
a in the isospin amplitudes are thought to be small for general reasons, the real 
parts of the b - coefficients (after taking a to be real) are determined by the well 
measured linear slopes g. However, it is easy to see from the isospin relations in 
tab.l that neither can the individual contributions bjj be over constraint (we 
have 3 measured linear slope parameters g for 3 quantities) nor is it possible 
to get any information on their imaginary parts or their phases, which are of 
primary interest here. In order to find these, we must take into account the 
quadratic slopes h, k. They are much more problematic from the experimental 
( statistics!) as well as from the theoretical (radiative corrections!) point of 
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view. 

After defining the phases /3 by 



6=| 6 | e if3 

it is easy to derive for each of the four channels a relation of the form 

h + U-?=4*pg. = R. (4) 
4cos 2 j3 | a | 2 

From the isospin components ajj, cjj of a, c given in table 1 it is clear that the 
r.h.s. R of this equation should be the same for the two charged channels and 
the two neutral channels respectively, independent of any assumptions on the 
phases of ajj, cij. It should be approximately the same for all four channels if 
C13 <C en , according to the AI = ^ rule. 

Experimental values for g, h, k, and R, calculated from (|j) assuming [3 = 0, 
are given in table 2. Instead of using PDG values, we choose to take here in each 
case the most significant (which in most cases means the latest) experiment. 
We would have to do this anyway for the 7r°7r°7r + and the 3 ir° data (which 
have not yet been included completely in the PDG tables). Moreover, because 
the consistency between some of the experimental data is certainly questionable 
(compare Q and Q), the errors may be more consistent; furthermore it becomes 
possible to compare our estimate for the effect of the Coulomb corrections with 
that given in the original paper Q. In any case, the PDG values do not disagree 
significantly from those used here. 

While the R - values from the two decay experiments are in perfect 
agreement, the situation for the charged K - decays is less clear, due to the 
lower (by one order of magnitude) statistics in the 7r°7r°7r + channel. If we assume 
c i3 *C en and take the R value from the (most significant) 7r 7r°7r° - experiment 
Q, we find linear relationships between h and k for all the other channels, as 
shown in fig.l. The measured (h, k) point should fall above the line, if (3 ^ 0. 
The consistency between the ir^n^n^ and the data is rather bad; one 

has to assume a fluctuation of ~ 2 st.dev., or a ratio of C13 / c\\ = 0.3 ± 0.1, 
in order to reach agreement. However from a comparison of figs. la, b which 
show the results from the same experiment derived without and with Coulomb 
corrections, respectively, it is clear that this result is very sensitive to these 
corrections. Therefore we will introduce more complete radiative corrections in 
the next section. 



Radiative Corrections 

Radiative corrections for hadronic processes are, generally speaking, model de- 
pendent, insofar as structural effects are concerned. The leading contributions 
can be estimated, however, by calculating the virtual photon exchanges and 
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soft photon emissions by point-like mesons, describing the weak decay by the 
most simple local interaction: 



L w = C W K 7Tl 7T2 7T3 . 



(5) 



To first order in a, the relevant graphs are shown in fig. 3a, while structural 
contributions fig. 3b are not considered. 

All analytic expressions needed are given below for easy reference. A FOR- 
TRAN code for calculating the correction factor as a function of X and Y for 
all K — > 3ir channels is available on request (schaale@ifh.de). 

The corrected decay probability is given by 



dT{ Si ) = dT ( Si )\l + — 



( a r 3 3 1 1 

L i=l i<j i=l J ^ 



(6) 



where ej = 0, ±1 are charge factors of K, irf, 



, [2Ae\ 2 f 1 , 
Zn 2 + -/n 

V fi J \ Vi 1 + Vi 



-i u-5 



— in 

2 Vl + ^i 



+ 2 



- 1 



l 



In 



4/i : 



1 - Qi(z) 



l + Qi(z) 



1-*) ln™-6^- ln™ + 
Si J n m — n fi 



and 
F 



In 



f2Ae\ 2 



+ (1 + 4) 



2 + 



+ 1 dz 



1 + % , 1 
fn- 



i ^ 2 -4 



in 



l + u t 



, 1 - Vi 
Uj,- tn : 

3 1 + Vi 

z 2 — v 2 



+ 2 



Z 2 V 2 j - 1 



—Qij(z) In 



1 - <2y(z) 



1 + 



where G = 0.915966 is Catalan's constant; 



1 , l-Vi 
— in — 2 ; 



f i 1 + Vi 

^Jz) = z 2 - 2z— - ( 1 - — ], T l = m 2 + n 2 , A = m 2 -fi 2 

Si \ Si J 



+ "" JW |(1 + + (1-^)^-1 

"^i = is the CMS-velocity; = y/l — 4/i 2 /sfc are the velocities in (i,j)- 

rest system; m = m K + , u = m w + , and Ae is the 7 cut-off energy. The conven- 
tional Coulomb correction factor || is given by 



Y[ Vij/ (exp(r]ij) - 1) , rjij = 2-Traeiej/\vi - Vj 

i<j 



(7) 



8G , 
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In fig. 4 are shown the correction factors according to (^), with Ae = WMeV 
are shown as functions of X, Y for all charged channels. They are largest for the 
— > decay (for this channel we show for comparison also the values 

calculated with Ae = 50MeV and with the conventional Coulomb factor ([?])). 
For the other decay modes the corrections are much less important and depend 
only on Y . For — ► 7r ± 7r°7r° the Coulomb factor (0) equals 1, because there 
is only one charged particle in the final state. 

We add some remarks concerning the well known singularities appearing 
in the treatment of radiative corrections, in order to show the limits of this 
approach. 

The first kind are ultraviolet divergences of the loop integrals which result 
in explicit cut-off dependent terms ~ ln(A/fi). The reason is that we have 
chosen purely local effective interactions without hadronic form factors, which 
would regularize these integrals. In our approach we renormalized them by the 
requirement that their contributions should disappear for transitions, in which 
the incoming charge reemerges as outgoing with the same velocity. This means 
we subtract FKn(si = (m — fi) 2 ) and F 1T1T {si = 0) respectively (standard renor- 
malization on mass shell)Q. As we neglected higher order corrections to L w 
(L = L w + rjLi, with r/ ~ 0.1, say), being responsible for the kinematic struc- 
ture of the Dalitz-plot, our results for slope parameters are correct up to terms 
of order arj. The constant terms au are affected by the unknown renormaliza- 
tion ambiguity of order a; in other words, hadron structure effects may differ 
between — > 3tt° and charged K decays by terms of order a. It is easy to 
show that this induces also corrections of order arj for slope parameters. Their 
treatment would require a more complicated effective Lagrangian, inclusion of 
structural photons etc., and would be strongly model-dependent. 

A second kind are so-called collinear singularities, appearing in the case of at 
least two charged particles with equal velocities in the final state, that means on 
the Dalitz-plot boundary. In these points the perturbation expansion breaks 
down. Accepting that some regions of the Dalitz-plot are just not handled 
correctly by the theory, some caution in the treatment of experimental data 
near these singularities may be required. It can be shown that no infinities are 
encountered in integrals over kinematic regions (integrability) . 

The third and last type are infrared divergences caused by low energy pho- 
tons. Their cure is well known, leading to the introduction of the upper limit 
Ae, defined here in the K rest system, up to which weak photons are to be 
included in the definition of the decay channel. A rough estimate of Ae may 
be derived from the mass resolution achieved in reconstructing the K mass. 
However, care should be exercised when high energy K-decays are analyzed. 
A correct experimental treatment would have to include radiative processes in 

lr The subtraction constants include also all contributions from photon emission graphs; in 
this respect the present results differ from those given in 113] 
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the Monte-Carlo and to establish an effective Ae in this way. 

In order to estimate the influence of the radiative and the Coulomb correc- 
tions on the Dalitz-plot parameters without reference to the experimental data, 
i.e. uncorrected Dalitz-plot densities, we introduce moments < X m ■ Y n >0 
with respect to normalized Dalitz-plot densities p(X,Y) and p'(X,Y), where 
p represents a constant density and p' includes radiative corrections to this 
constant density. After expanding p' in X, Y, it is easy to get approximate 
expressions for p' moments < ... >' in terms of p moments < ... >, but for our 
case we have calculated all relevant moments numerically. If we suppose that 
the experimental data sample is already corrected for experimental efficiency 
and background []. , the slope parameters are to be derived from a fit to the 
Dalitz-plot density, which may be written as 

where for convenience the vectors 

/ Y 

v= Y 2 
\ X 2 

are introduced. We write v(i) = v{Xi,Y.{) for measured Xi and Yi for the 
i th event, g is to be estimated by the Maximum Likelihood method from the 
Likelihood function 

n 

L = Hf(Xi,Yi) (9) 

i=i 

leading to the system of equations for g 

iv % l*^ (io) 

n t-r 1 1 + v(i) ■ g 1+ < v > g 

(we do not distinguish here between estimates and population values). If we 
now apply a further correction, e.g. p — > p' , we get, with the same sample of 
experimental data Xi, Y^, corrected parameters g',h',k', which can be expressed 
in terms of g and moments < ... >, < ... >'. The equations for <f are: 

1 v(i) <v>' 




+ v(i) ■ g' 1+ < v >' -g' 



2 There are significant differences, at least concerning the parameter k, between |6j and a 
later experiment II, which found k = -0.0205 ±0.0039 to be compared with -0.0075 ±0.0019 
|j] for K + — » 7r + 7T + 7r~. As the first experiment measured only the momentum of the odd 
pion in the final state, the second one all momenta, effects of the above mentioned kind may 
be present. |)j do not present data without (Coulomb) correction, therefore we used only 

3 In the actual evaluation of an experiment, p would be taken from a Monte-Carlo simulation 
of the measured Dalitz-plot distribution, using a constant Dalitz-plot density as input. 
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For small corrections we may expand both sides in Ag = g 1 — g and obtain a 
linearized set of equations: 

V ■ Ag = <j> „- <*_>', „ (12) 
l+<v>-g l+<v>'-g y ' 

where the symmetric matrix V has the elements 

\r 1 \ - ( VkVl \ < Vk >'< vi >' , , Q s 

Vki = - }^ [ n | ^ 2 L ~ n , . -, , , -» 2 [k, 1 = 1,2,3). (13) 
n V (1 + v ■ g) £ > v=v(%) (1+ < v > -gY 

The dependence on X{, can now be eliminated by replacing ^ J2 with / dXdY f(X, Y), 
leading to 

V 1 / V JM \ < y k >'< v _i >' {U \ 

kl 1+ <v> -g\l + v-g/ (l+<v>'-g) 2 ' { 1 

where the first terms on the right-hand side are also to be evaluated for given 
g,h,k numerically. 

Surely one has to be aware of the severe limitations of the above approach if, 
for experimental or other reasons, the uncorrected parameters g do not represent 
the density over the whole Dalitz-plot. To give an extreme example, suppose 
they had been derived from a fit to the Dalitz-plot density in a region where 
the corrections disappear. They were then found identical to the corrected 
parameters, if the "corrections" are applied to the raw data sample. For our 
method one has to assume however that the uncorrected parameters fit the 
density equally well over the whole Dalitz-plot. If, as we may further suppose 
in our example, there are sizeable corrections to the density in the unmeasured 
region this is not the case, and consequently we find the corrected parameters 
different from the uncorrected ones, possibly even outside the statistical errors. 

As a check of our method we calculated for the conventional Coulomb cor- 
rection factor the corrected parameters if, corresponding to the first column of 
table 2, from the uncorrected values of ref.|6| in the second column. The results 
are given in table 3 together with the differences with respect to the corrected 
values from ref . Q . 

Our conclusion from this comparison is that the method is useful to demon- 
strate the influence of radiative corrections on the quadratic slope parameters, 
where statistical errors are relatively large. It is not a substitute, however, for 
a complete (re) analysis of precision experimental data like those existing for 
linear slope parameters, for which systematic corrections are more subtle. 



Results and Conclusions 

The results for Ag, Ah, Ak are shown as a function of the linear slope g in 
fig. 5 and for the actual g-values (see table 2) in table 4. The dependence on 
h and k is for small values of these parameters negligible, they are set to zero 
everywhere. 
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For the if ± — ► ir :il TT + TT channel the discrepancy of the i?-value with that 
for K® —>■ 37T° (see table 2) disappears after applying the radiative corrections 
(||) instead of the Coulomb-factor (^(compare fig. la) with fig. 2). We find 

^xo,r.(/3 = 0) = --0078 ±.0089, 

indicating that the quadratic coefficients c in the decay amplitude (||) may be 
small and of the same order as for K\ -decays. 

The corresponding corrections for the channel — > n^n^ir® are smaller 
than the numerical accuracies of the numbers in table 2. A comparison of 
i?-values of the different charged K-decay channels with comparable statistics 
would be interesting. 

In order to derive some first information on possible phases (3 we may use 
the relation 

a 2 

R{(3) = R(0) - ^-tan 2 /?. 

For the K° L — > 7r + 7r 7r°-channel we can identify R + with R 000 and find 
from the last row of tab. 20 

tan 2 (3 = 0.0 ± .086 , 

i.e. a phase angle (3^ < 16°. For the charged kaon decay — ► 7r ± 7r + 7r _ , 
assuming C13 <C c\\ and comparing R 000 also with Rf^~ orr , we find analogously 
(3 ±+ - < 43°. 

Despite comparable errors of the i?-values for the two cases, the restriction 
on /? ±_l is weaker than that found on f3^ °. This is due to the different 
linear slopes g. In view of this, and also because it is much less influenced 
by radiative corrections, the channel — > 7r ± 7r 7r (with g ~ 0.6) deserves 
special attention by experiment. Enhancing the data sample for this channel 
by an order of magnitude (to ~ 5 • 10 5 events) could lead to a determination of 
the phase f3 with an error < 15° similar as for K® decays. Clearly, this would 
help a lot to constrain effective Lagrangian models with regard to higher order 
(p 4 -, loop-, penguin-) contributions, especially if taken together with results 
from radiative K-decays, where also new experimental and theoretical work is 
going on. 

Besides this, one should be aware that also for the other channels consid- 
ered here one has to rely presently on only few large statistics experiments. 
Concerning the presentation of new data, we would like to advocate to publish 
the data also in a form uncorrected for radiative or Coulomb effects. Further- 
more, one should clearly state the regions of the Dalitz-plot to which the slope 
parameters have been fitted. There is room for later improvements of the ra- 
diative corrections, including structural radiation, after more realistic effective 
Lagrangians will have been introduced. 

4 radiative corrections by the authors [hoi are already included 
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Channel 
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00+ 


+ 


- 


000 


a 


2(an + ai 3 ) 


an + 013 


-(an 


- 2ai 3 ) 


-3(on - 2ai 3 ) 


b 


-(6n + 613) + 623 


fen + fel3 + fe23 


-(fen 


- 26i 3 ) 





c 


2(cii + ci 3 ) 


Cll + C13 


-(en 


- 2ci 3 ) 


-3(cn - 2c 13 ) 


d 


-((in + dl 3 ) +^23 


dn + d 13 + d 23 


-(d n 


- 2d 13 ) 






Table 1: Isospin Amplitudes 



Ch. 




±7TF 


TT TT°Tr + 


7T + 7T 7T° 




Expt. 


§,a 




Q 




i 


9 


-.2173 ± .0026 


-.1866 + .0025 


.575 ± .022 


.677 ± .010 





h 


.0156 ± .0062 


.00125 ± .0062 


.021 ± .023 


.079 ± .007 


-.0033 ± .0013 


k 


-.0079 ± .0019 


.0029 ± .0021 


.011 ± .007 


.0097 ± .0018 


h/3 


R 


-.0199 ± .0084 


.0013 ± .0088 


-.029 ± .032 


-.0065 db .0095 


-.0066 ± .0026 



Table 2: Experimental Data (a with, b without Coulomb correction) 





calc. by (0) 


Diff. to § 




9 


-.2236 ± .0025 


.0063 ± .0037 




h 


.0149 ± .0062 


.0007 ± .0088 




k 


-.0079 ± .0021 


0.0 ± .003 



Table 3: Comparison for Coulomb corrections 



Channel 


± + - 


±00 


+ -0 


Ao 


-1.99- 10" 2 


-1.45 ■ 10" 3 


7.02 • 10" 3 


Ah 


7.28 • 10~ 3 


-8.47 • 10" 4 


1.63- 10~ 3 


Ak 


-4.72 • 10" 3 


7.66 • 10~ 7 


-1.70 • 10" 5 



Table 4: Radiative Corrections for Slope Parameters 
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Figure captions 



Fig. 1 Plots of quadratic slope parameters k vs. h from table 2. The linear 
relation (||) is indicated with 1 s.d. errors. 
+ H — a : channel — > Tr^vr+Tr - with Coulomb corrections (|?j) 
uncorrected 

channel -fT 1 * 1 — ► vr^vr ^ , uncorrected 



b 

+ 00 
+ - 



channel — ► 7r + 7r 7r°, corrected by [10] 



Fig. 2 The same as fig.l for channel — > it^^n corrected according to (|T^) 
with rad. corr. (M) 



Fig. 3 a) Graphs for radiative corrections to first order 
b) Graphs with inclusion of structural radiation 

Fig. 4 Plots of correction factors as functions of X, Y 

a) Rad. corr. for -> ir^ir+ir' , Ae = lOMeV, 50MeV (broken lines) 

b) Coul. corr. for K — * tt^tt^tt^ 

c) Rad. corr. for — ► ir + ir~ir° 

d) Rad. corr. for — > it^it^-K® 

Fig. 5 a) - c) Dependence of Ag, Ah, Ak on g (for h = k = 0) for the channel 
— y 7r ± 7r + 7r~ 



13 



This figure "figl-l.png" is available in "png" format from: 
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This figure "fig2-l.png" is available in "png" format from: 
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